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1.  INTRODUCTION 

In  recent  years  composites,  especially  f iber~reinforced 
laminates,  have  found  increasing  application  in  many 
engineering  structures.  This  is  mainly  due  to  two  desirable 
features  of  fiber-reinforced  composites:  (i)  the  high 
stiffness-to-weight  ratio,  and  (ii)  the  anisotropic  material 
property  that  can  be  tailored  through  variation  of  the  fiber 
orientation  and  stacking  sequence  of  lamina  -  a  feature  that 
gives  the  designer  an  added  degree  of  flexibility.  The 
increased  use  of  fiber-reinforced  composites  as  structural 
elements  has  generated  considerable  interest  in  the  analysis 
of  lamianted  (anisotropic)  composite  plates. 

Recent  developments  in  the  analysis  of  plates  laminated 
of  fiber-reinforced  materials  indicate  that  the  plate 
thickness  has  more  pronounced  effect  on  the  behavior  of 
composite  plates  than  isotropic  plates.  The  classical  thin- 
plate  theory  (CPT)  assumes  that  normals  to  the  midsurface 
before  deformation  remain  straight  and  normal  to  the 
midsurface  after  deformation,  implying  that  thickness  shear 
deformation  effects  are  negligible.  As  a  result,  the  natural 
frequencies,  for  example,  calculated  using  the  thin-plate 
theory  are  higher  than  those  obtained  by  including  the 
transverse  shear  deformation  effects.  Also,  the  transverse 
deflections  predicted  by  the  thin-plate  theory  are  lower  than 
those  predicted  by  a  shear  deformable  theory  (SOT).  Due  to 
the  low  transverse  shear  modulus  relative  to  the  in-plane 
Young's  moduli,  the  transverse  shear  deformation  effects  are 
even  more  pronounced  in  the  composite  plates.  Thus  a  reliable 
prediction  of  the  small  deflection  response  characteristics  of 
high  modulus  composite  plates  requires  the  use  of  shear 
deformable  theories. 
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When  the  transverse  deflections  experienced  by  plates  are 
not  small  in  comparison  with  the  plate  thickness,  the 
interaction  between  membane  stresses  and  the  (bending  and 
shear)  curvatures  of  the  plate  must  be  considered.  The 
interaction  results  in  midplane  stretching,  which  leads  to 
nonlinear  terms  in  the  equations  of  motion.  Thus,  a  more 
acccurate  prediction  of  deflections,  stresses,  and  frequencies 
requires  the  solution  of  the  plate  equations  that  account  for 
large  deflections  and  thickness  shear  deformation. 

The  present  paper  is  concerned  with  a  review  of  recent 
developments  in  the  finite-element  analysis  of  layered 
composite  plates.  The  paper  presents  a  shear  deformable 
theory  that  accounts  for  large  rotations  (in  the  von  Karman 
sense),  a  finite  element  model  of  the  equations,  and  finite 
element  results  (of  linear  and  nonlinear  theory)  for  static 
bending  and  free  vibration  of  rectangular  plates.  The 
following  brief  review,  which  cites  the  papers  that  have 
appeared  in  the  open  literature  on  the  linear  and  nonlinear 
bending  and  vibration  of  layered  composite  plates,  provides  a 
background  for  the  present  paper. 

A  number  of  shear  deformable  theories  for  laminated 
plates  have  been  proposed  to  date.  The  first  such  theory  for 
laminated  isotropic  plates  is  due  to  Stavsky  [1].  The  theory 
has  been  generalized  to  laminated  anisotropic  plates  by  Yang, 
Norris,  and  Stavsky  [2].  The  Yang-Norris-Stavsky  (YNS)  theory 
represents  a  generalization  of  Reissner-Mindl in  plate  theory 
for  homogeneous,  isotropic  plates  to  arbitrarily  laminated 
anisotropic  plates  and  includes  shear  deformation  and  rotatory 
inertia  effects.  A  review  of  various  other  theories,  for 
example,  the  effective  stiffness  theory  of  Sun  and  Whitney 
[3],  the  higher-order  theory  of  Whitney  and  Sun  [4],  and  the 
three-dimensional  elasticity  theory  of  Srinivas  et  al.  [5-7], 
can  be  found  in  the  paper  by  Bert  [8].  It  has  been  shown 
(see,  for  example,  [3,5,9-13])  that  the  YNS  theory  is  adequate 
for  predicting  the  overall  behavior  such  as  transverse 
deflections  and  natural  frequencies  (first  few  modes)  of 
laminated  anisotropic  plates. 

The  first  application  of  the  YNS  theory  is  apparently  due 
to  Whitney  and  Pagano  [14],  who  presented  closed-form 
solutions  for  symmetric  and  antisymmetric  cross-ply  and  angle- 
ply  rectangular  plates  under  sinusoidal  load  distribution  and 
for  free  vibration  of  antisymmetric  angle-ply  rectangular 
plates.  Fortier  and  Rosettos  [15]  analyzed  free  vibration  of 
thick  rectangular  plates  of  unsymmetric  cross-ply  construction 
while  Sinha  and  Rath  [16]  considered  both  vibration  and 
buckling  for  the  same  type  of  plates.  Following  Whitney  and 
Pagano  [14],  Bert  and  Chen  [17]  presented  closed-form  solution 
for  the  free  vibration  of  simply  supported  rectangular  plates 
of  antisymmetric  angle-ply  laminates. 


Finite-element  analysis  of  layered  composite  plates  began 
with  the  works  of  Pryor  and  Barker  [18],  and  Barker,  Lin  and 
Dana  [19],  who  employed  an  element  with  seven  degrees  of 
freedom  (three  displacements,  two  rotations  and  two  shear 
slopes)  per  node  to  analyze  thick  laminated  plates.  Mau, 

Tong,  and  Pi  an  [20],  and  Mau  and  Witmer  [21]  employed  the  so- 
called  hybrid-stress  finite  element  method  to  analyze  thick 
composite  plates.  Using  finite  element  models  based  on  a  form 
of  Reissner's  plate  theory  (i.e.,  'mixed'  formulation),  Noor 
and  Mathers  [22,23]  studied  the  effects  of  shear  deformation 
and  anisotropy  on  the  response  of  lamiated  anisotropic  plates. 
One  of  the  elements  used  by  the  authors  (see  [22])  possesses 
eighty  degrees  of  freedom  per  element,  and  one  cannot  afford 
to  use  it  without  the  luxury  of  having  the  enormous 
computational  resources  the  authors  had.  Hinton  [24,25] 
employed  the  so-called  finite  strip  method  to  study  the 
flexure  and  free  vibration  of  layered  cross-ply  plates. 

Mawenya  and  Davies  [26],  and  Panda  and  Natara.ian  [27]  used  the 
quadratic  shell  element  of  Ahmad,  Irons  and  Zienkiewicz  [28] 
to  analyze  bending  of  thick  multilayer  plates.  Using  two 
hybrid-stress  elements,  Spilker,  Chou,  and  Orringer  [29] 
studied  the  static  bending  of  layered  composite  plates.  The 
number  of  degrees  of  freedom  in  one  of  the  elements  is 
proportional  to  the  number  of  layers,  and  therefore  the  core 
storage  and  execution  time  requirements  for  the  element 
increase  rapidly  with  the  number  of  layers  in  the  plate. 
Recently,  the  present  author  developed  a  simple  and  efficient 
finite  element  based  on  the  YNS  theory  [30].  The  element 
contains  three  displacements  and  two  shear  slopes  as  degrees 
of  freedom  per  node.  The  accuracy  and  convergence 
characteristics  of  the  element  were  investigated  in  [31].  The 
element  has  been  used  successfully  in  the  thermoelastic 
analysis  of  ordinary  and  bimodulus  (i.e.,  different  elastic 
properties  in  tension  and  compression)  layered  composite 
plates  by  the  author  and  his  colleagues  [30-34]. 

Much  of  the  previous  research  in  the  analysis  of 
composite  plates  is  limited  to  linear  problems.  This  may  have 
been  due  to  the  complexity  of  the  nonlinear  partial 
differential  equation  associated  with  the  large-deflection 
theory  of  composite  plates.  Approximate  solutions  to  the 
large-deflection  theory  (in  von  karman's  sense)  of  laminated 
composite  plates  were  attempted  by  Whitney  and  Leissa  [35], 
Bennett  [36],  Bert  [37],  Chandra  and  Raju  [38,39],  Zaqhloul 
and  Kennedy  [40],  Chia  and  Prabhakara  [41,42],  and  Noor  and 
Hartley  [43].  Chandra  and  Raju  [38,39],  and  Chia  and 
Prabhakara  [41,42]  employed  the  Galerkin  method  to  reduce  the 
governing  nonlinear  partial  differential  equations  to  an 
ordinary  differential  equation  in  time  for  the  mode  shape;  the 
perturbation  technique  was  used  to  solve  the  resulting 
equation.  Zaqhloul  and  Kennedy  [40]  used  a  finite-difference 
successive  iterative  technique  in  their  analysis.  In  all  of 
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these  studies,  with  the  exception  of  [43],  the  transverse 
shear  effects  were  neglected.  The  finite  element  employed  by 
Noor  and  Hartley  [43]  includes  the  effect  of  transverse  shear 
strains;  however,  it  is  algebraically  complex  and  involves  a 
large  degree  of  freedom  per  element  and  thus  one  can  preclude 
the  use  of  such  elements  in  the  nonlinear  analysis  of 
composite  plates.  Recently,  the  finite  element  developed  in 
[30]  by  the  author  was  extended  to  nonlinear  bending  of 
composite  plates  [44,45]. 

Analysis  of  nonlinear  vibration  of  single-layer 
orthotropic  plates  was  considered  by  Ambartsumyan  [46],  and 
Hassert  and  Nowinski  [47].  Nowinski  [48,49]  analyzed 
recti  1 inearly  orthotropic  plates  of  circular  and  triangular 
planforms  using  Galerkin  method.  In  these  studies  the  effect 
of  transverse  shear  deformation  and  rotatory  inertia  were  not 
considered.  Wu  and  Vinson  [50]  presented  the  dynamic  analogue 
of  Berger's  equation  of  motion  for  an  orthotropic  plate, 
including  the  effect  of  transverse  shear  deformation  and 
rotatory  inetia;  however,  the  solutions  were  restricted  to 
only  the  transverse  shear  deformation.  Mayberry  and  Bert  [51] 
presented  experimental  as  well  as  theoretical  work  on 
nonlinear  vibration  of  laminated  plates;  however,  the 
theoretical  investigation  was  limited  to  single-layer 
specially  orthotropic  rectangular  plate  with  all  four  edges 
clamped,  and  not  including  the  effect  of  transverse  shear 
deformation  and  rotatory  inertia.  Using  an  assumed  mode  shape 
and  Galerkin  method,  Nowinski  [52]  presented  a  general 
equation  for  the  nonlinear  analysis  (i.e.,  large  deflection 
and  large  amplitude  free  vibration)  of  orthotropic  plates  (see 
also,  Sathyamoorthy  and  Pandalai  [53]).  Prabhakara  and  Chia 
[54]  presented  an  analytical  investigation  of  the  nonlinear 
vibration  of  a  rectangular  orthotropic  plate  with  all  simply 
supported  and  all  clamped  edges.  The  effect  of  transverse 
shear  and  rotatory  inertia  on  large  amplitude  vibration  of 
composite  plates  was  reoorted  recently  by  Sathyamoorthy  and 
Chia  [55-57]  using  Galerkin  method  and  the  Runge-Kutta 
numerical  procedure. 

Layered  composite  plates  exhibit,  in  general,  coupling 
between  the  in-plane  displacements  and  the  transverse 
displacement  and  shear  rotations.  For  plates  having  layers 
stacked  symmetrically  with  respect  to  the  mid-plane,  the 
bending-stretching  coupling  terms  vanish  and  the  problem 
becomes  relatively  simpler.  Wu  and  Vinson  [58]  extended  their 
earlier  work  [50]  to  deal  with  the  nonlinear  vibration  of 
symmetrically  stacked  laminated  composite  plates.  The  first 
nonlinear  vibration  analysis  of  unsymmetrical ly  laminated 
plates  is  due  to  Bennett  [36],  who  considered  simply  supported 
(with  immovable  edges)  angle-ply  plates.  Using  the  thin  plate 
theory  of  layered  composite  plates  and  Galerkin  method,  Bert 
[37]  investigated  the  nonlinear  vibration  of  a  rectangular 
plate  arbitrarily  laminated  of  anisotropic  material.  A 
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multimode  (two-term)  solution  for  nonlinear  vibration  of 
unsymmetric  all-clamped  and  all-simply  supported  angle-ply  and 
cross-ply  laminated  plates  was  reported  by  Chandra  and  Basava 
Raju  [38,39,59].  Chandra  [39]  used  a  one-term  Galerkin 
approximation  for  the  dynamic  von  Karman  plate  equations  and 
the  perturbation  technique  for  the  resulting  ordinary  equation 
in  time  to  investigate  the  large-amplitude  vibration  of  a 
cross-ply  plate  which  is  simply  supported  at  two  opposite 
edges  and  clamped  at  the  other  two  edges.  Prabhakara  and  Chi  a 
[54]  presented  an  analytical  investigation  of  the  nonlinear 
free  flexural  vibrations  of  unsymmetric  cross-ply  and  angle- 
ply  plates  with  all -clamped  and  all-simply  supported  edges. 

The  normal  and  tangential  boundary  forces  in  the  plane  of  the 
plate  were  assumed  to  be  zero,  ^n  all  of  these  studies,  the 
effects  of  shear  deformation  and  rotatory  inertia  were 
neglected,  and  only  rectangular  geometries  were  considered. 

The  latter  restriction  is  a  direct  result  of  the  methods  of 
analysis  used  (i.e.,  Galerkin  method,  perturbation  method  and 
the  double  series  method  cannot  be  applied  to  plates  of 
complicated  geometries).  Recently,  the  present  author  [60] 
investigated  the  nonlinear  vibration  of  layered  composite 
plates  using  the  finite  element  method. 

2.  GOVERNING  EQUATIONS  ([45,60]) 

Consider  a  (thick)  plate  of  uniform  thickness  h, 
laminated  of  anisotropic  layers  with  the  material  axes  of  each 
layer  being  arbitrarily  oriented  with  respect  to  the  midplane 
of  the  plate.  The  midplane  of  the  plate,  denoted  by  R,  forms 
the  x-y  plane  of  the  coordinate  system,  with  the  z-axis 
perpendicular  to  the  midplane  of  the  plate.  In  order  to 
account  for  the  midplane  stretching  due  to  large  deflections, 
the  shear  deformable  theory  of  Whitney  and  Pagano  [14]  is 
modified  to  include  large  rotations  (in  the  von  Karman  sense). 
The  displacement  field  is  assumed  to  be  of  the  form, 

Uj(x,y,z,t)  =  u(x,y,t)  +  z  4'x(x,y,t), 

U2(x,y,z,t)  =  v(x,y,t)  +  z  4'y(x,y,t),  (1) 

U3(x,y,z,t)  =  w(x,y,t). 

Here  t  is  the  time;  u^,  U2.  Uj  are  the  displacements  in  x,  y, 
z  directions,  respectively;  u,  v,  w  are  the  associated 
midplane  displacements;  and  and  i|;y  are  the  slopes  in 
the  xz  and  yz  planes  due  to  bending  only.  Assuming  that  the 
plate  is  moderately  thick  and  strains  are  much  smaller  than 
rotations,  the  nonlinear  strain-displacement  relations  can  be 
expressed  in  the  for.i, 
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wherein  the  products  of  iii^,  liiy,  3Uj/3x  and  3u2/3y  are 
neglected.  Since  the  constitutive  relations,  to  be  given 
shortly,  are  based  on  the  plane-stress  assumption,  strain  Ej 
does  not  come  into  the  eguations. 


Neglecting  the  body  moments  and  surface  shearing  forces, 
we  write  the  equations  of  motion  (in  the  absence  of  surface 
and  body  forces)  as. 
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where  P,  R,  and  I  are  the  normal,  coupled  normal -rotatory,  and 
rotatory  inertia  coefficients. 
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p(m)  being  the  material  density  of  the  m-th  layer,  N^,  Q^,  and 
are  the  stress  and  moment  resultants  defined  by 
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and  N(*)  is  the  nonlinear  operator, 

N(w,Ni )  =  (Nj  -^)  +  -^  (N^  -^)  +  lx  (’'*6 

('^2  -37) 

Here  (i  =  1,2... 6)  denotes  the  in-plane  stress  components 
■  ®X’  °2  "  °y’  ‘’■i  ”  °yz’  ■  ®xz  °fa  “  ‘’xy)  • 

Assuming  the  existence  of  one  plane  of  elastic  symmetry 
parallel  to  the  plane  of  each  layer,  the  constitutive 
equations  for  the  m-th  layer  (in  the  plate  coordinates)  are 
given  by 


[Q 


where  are  the  stiffness  coefficients  of  the  m-th  layer  in 
the  plate  coordinates.  Combining  eqns.  (5)  and  (6),  we  obtain 
the  plate  constitutive  equations, 

lOj)  LA45A5j^J|£5| 


The  A^j,  B^j,  D^j  (i,j  =  1,2,6),  and  Ajj  (i,j  =  4,5)  are  the 
respective  inplane,  bending-plane  coupling,  bending  or 
twisting,  and  thickness-shear  stiffnesses,  respectively: 

^Ai j’^i j’^i j)  ■  ^  Qi j^(l.z.2^)dz  , 

*'  ^m 

Here  denotes  the  distance  from  the  mid-plane  to  the  lower 
surface  of  the  m-th  layer,  and  are  the  shear  correction 

coefficients. 

Equations  (3)  and  (7)  must  be  adjoined  by  appropriate 
boundary  conditions  of  the  problem.  The  variational 
formulation  of  these  equations  indicate  the  following 
essential  and  natural  boundary  conditions; 

essential:  specify,  Un,  u^,  w,  i|(„,  ij;,, 

n  5  n  s 

natural;  specify,  N^,  q, 

wherein  Up  and  Uj,  for  example,  denote  the  normal  and 
tangential  components  of  the  in-plane  displacement  vector,  u,  = 
(u,v). 

3.  FINITE  ELEMENT  MODEL  ([44,45,60]) 

The  equations  presented  in  the  previous  section  are  valid 
for  any  arbitrarily  laminated  composite  plate.  Except  under 
certain  special  conditions  of  geometry,  boundary  conditions 
and  stacking  of  layers,  these  equations  do  not  admit  exact 
solutions  (in  fact,  exact  closed-form  solutions  cannot  be 
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obtained  without  invoking  smal 1 -displacement  assumption;  see 
[29]).  Here  we  present  a  finite-element  model  associated  with 
the  equations  governing  the  nonlinear  analysis  of  layered 
composite  plates.  The  element  is  an  extension  of  the  penalty 
plate-bending  element  developed  for  the  linear  analysis  of 
layered  composite  plates  by  the  author  [28,29]. 

Consider  a  finite-element  analogue,  IT,  of  the  midplane  of 
the  plate,  R.  Over  a  typical  element  Rg  of  the  mesh  IT  a 
typical  generalized  displacement,  U,  is  interpolated  in  space 
by  expression  of  the  form, 

U  =  2  (10) 

i  '  ' 


where  Uj  is  the  value  of  U  at  node  i,  ij>j  is  the  finite-element 
interpolation  function  at  node  i  and  n  is  the  number  of  nodes 
in  the  element.  For  simplicity,  we  use  the  same  interpolation 
for  each  of  the  generalized  displacements. 


The  weak  form  of  eqns.  (3)  and  (7)  over  a  typical  element 
is  given  by 

0  =J^e^^u{Pu^tt+R'^'x,tt)'■^^xNl+«u,yN6+«v(Pv^U+R'^'y,tt) 

+  6v^^N6+6v^yN^+6w  Pw^tt)^^«,xQ]+^,yQ2 

a.  36w  3w  M  j.36w  3w  M  +  36w  3w  u  g. 36w  3w  u 
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+  «'i'y(I'l'y,tt'^Rv^tt)-*'^'*'y,xM6+^'l'y,/2+^'<'yQ2^ 

*  f.  («ll„N„«U5»„5)dst/  6w,ds+/  (6»„M„«*5H„5)ds, 

(11) 


where  N^,  M^,  and  are  given  in  terms  of  the  generalized 
displacements  by  eqn.  (7).  The  boundary  terms  N^,  q,  Mp, 
and  Mpj  (notation  used  is  very  standard)  in  eqn.  (ll)  get 
cancelled  at  interelement  boundaries,  and  equal  to  those 
specified  at  the  plate  boundary. 

Let  u,  V,  w,  and  be  approximated  over  each  element 

t>y  ^ 

U  =  UT(t),  V  =  VT(t),  W  =  WX(t),  l|;^  =  Xu(t),  =  Yu(t) 

(12) 

where  U,  V,  etc.,  are  given  by  ean,  (10),  and  r(t),  x(t),  and 
M(t)  are  time  dependent  functions  whose  specific  form  is  to 


be  determined.  Substituting  eqns.  (10)  and  (12)  into  eqn. 
(11),  we  obtain  the  following  element  equation: 


[K]{A}  +  [M]{A}  =  {F}  (13) 

Here  {A}  is  the  column  vector  of  the  nodal  values  of  the 
generalized  displacements,  [K]  is  the  matrix  of  stiffness 
coefficients  (which  depends  on  the  generalized  displacements), 
[M]  is  the  matrix  of  mass  coefficients,  and  (F)  is  the  column 
vector  containing  the  boundary  contributions.  The  elements  of 
[K]  and  [M]  are  given  in  Appendix  B.  It  should  be  noted  from 
[K]  in  the  appendix  that  the  time  functions  are  not  harmonic. 
That  is,  strictly  speaking,  eqn.  (13)  must  be  soT ved  as  one 
of  transient  equation  (even  in  the  case  of  free  vibrations). 
However,  in  the  present  analysis  we  assume,  for  simplicity, 
that 


T  =  u  =  =  COS^  u)t 


(14) 


and  retain  only  the  first  term  of  the  cosine  series.  This 
assumption  yields  the  standard  eigenvalue  problem  in  the  case 
of  natural  vibration: 


([K]  -  w2[m]){a}  =  {0}. 


(15) 


The  solution  procedure  consists  of  a  direct  iteration,  in 
which  the  global  stiffness  [K]  is  updated  using  the  global 
displacement  (eigen  function)  vector  (A)  from  the  previous 
iteration.  If  {A}  is  set  to  zero  (as  was  done  in  the  present 
analysis)  at  the  beginning  of  the  iteration  procedure,  we 
obtain  the  linear  solution  (frequencies)  of  the  problem  at  the 
end  of  the  first  iteration.  The  iteration  is  terminated  when 
the  nonlinear  solution  (frequencies)  obtained  during  two 
consecutive  iterations  differ  by  some  small  number  (say, 
one  percent). 

In  the  present  study  linear  (4-node)  and  quadratic  (8- 
node  and  9-node)  isoparametric  elements  are  employed.  The 
element  stiffness  matrices  are  of  the  order  20  x  20,  40  x  40, 
and  45  x  45  respectively. 


As  pointed  out  in  [30,61],  the  shear  deformable  theory 
presented  herein  can  be  derived  from  the  classical  thin  plate 
theory  by  using  the  penalty  function  method  to  incorporate  the 
slope-displacement  relations 


3w  _ 
■ 


as  constraints  into  the  variational  formulation  of  the  thin 
plate  theory.  H  is  well-known  that  reduced  integration  must 
be  used  in  the  numerical  solution  of  penalty  function 
problems.  In  other  words,  the  shear  energy  terms  (which 
correspond  to  the  penalty  functional)  in  the  element  matrices 
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must  be  underinteqrated.  This  means  that  a  1x1  Gauss  rule 
must  be  used  in  place  of  the  standard  2x2  Gauss  rule  in  the 
evaluation  of  the  stiffness  coefficients  associated  with  the 
shear  energy  for  the  four-node  element.  For  more  details  on 
the  reduced  integration  in  plates,  see  [62,63]. 

6.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  finite  element  presented  herein  was  employed  in  the 
nonlinear  analysis  of  rectangular  plates.  The  following 
material  properties  typical  of  advanced  fiber-reinforced 
composites  were  used  in  the  present  study: 

Material  I;  E^/E2=25,  G^2/E2=0-5,  G22/E  =0.2,  Vj2~0-25 

Material  II;  E^/E2=40,  612/^2=  0.6,  623/^2=  0.5, 

(16) 

It  is  assumed  that  6^3  =  G23  and  v^2  ■  ^13*  4  value  of  5/6 

was  used  for  the  shear  correction  coefficients, 

Whitney  [64]).  All  of  the  computations  were  carried  on  an  IBM 
370/158  computer  in  double  precision. 

6.1  Static  Bending  Results  ([30,31,44,45]). 

To  show  the  effect  of  the  reduced  integration,  and  to 
illustrate  the  accuracy  of  the  present  element,  results  of  the 
linear  analysis  are  presented.  Numerical  experiments  were 
conducted  to  investigate  the  effect  of  element  type  (i.e., 
linear,  eight-node  quadratic,  and  nine-node  quadratic),  mesh 
(L2  =  linear  2  by  2,  Q2  =  2Q8  =  8-node  quadratic  2  by  2 
elements  in  quarter  plate),  and  reduced  integration  on  the 
deflections  and  stresses. 

In  Table  1  linear  and  quadratic  elements  are  compared  for 
deflections  and  stresses  of  a  three-layer  (hj  =  63  =  h/4,  h2  = 
h/2)  square  plate  (material  I)  subjected  to  sinusoidal 
loading.  This  problem  is  also  equivalent  to  a  four-layer 
(equal  thickness)  cross-ply  laminate.  The  following  boundary 
conditions  were  used  (in  the  finite  element  method  only  the 
essential  boundary  conditions  were  imposed): 

Uo(x,0)  =  Uo(x,b)  =  0,  N2(x,0)  =  N2(x,b)  =  0 

Vo(0*y)  =  V(,(a,y)  =  0,  Nj(0,y)  =  Ni(a,y)  =  0 

SS-1:  w{x,0)  =  w(x,b)  =  w(0,y)  =  w(a,y)  =  0 

if'x(x,0)  =  ij'x(x,b)  =  0,  M2(x,0)  =  M2(x,b)  =  0 

'l'y(0,y)  =  'l'y(a,y)  =  0,  Mj(n,y)  =  Mi(a,y)  =  0 


(17) 
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be  determined.  Substituting  eqns.  (10)  and  (12)  into  eqn. 
(11),  we  obtain  the  following  element  equation; 

[K]{A}  +  [M]{A}  =  {F}  (13) 

Here  {A}  is  the  column  vector  of  the  nodal  values  of  the 
generalized  displacements,  [K]  is  the  matrix  of  stiffness 
coefficients  (which  depends  on  the  generalized  displacements), 
[M]  is  the  matrix  of  mass  coefficients,  and  {F}  is  the  column 
vector  containing  the  boundary  contributions.  The  elements  of 
[K]  and  [M]  are  given  in  Appendix  B.  It  should  be  noted  from 
[K]  in  the  appendix  that  the  time  functions  are  not  harmonic. 
That  is,  strictly  speaking,  eqn.  (13)  must  be  solved  as  one 
of  transient  equation  (even  in  the  case  of  free  vibrations). 
However,  in  the  present  analysis  we  assume,  for  simplicity, 
that 


T  =  u  =  =  cos^  cot  (14) 

and  retain  only  the  first  term  of  the  cosine  series.  This 
assumption  yields  the  standard  eigenvalue  problem  in  the  case 
of  natural  vibration: 


([K]  -  oo2[M]){a}  =  (0).  (15) 

The  solution  procedure  consists  of  a  direct  iteration,  in 
which  the  global  stiffness  [K]  is  updated  using  the  global 
displacement  (eigen  function)  vector  (a)  from  the  previous 
iteration.  If  {A}  is  set  to  zero  (as  was  done  in  the  present 
analysis)  at  the  beginning  of  the  iteration  procedure,  we 
obtain  the  linear  solution  (frequencies)  of  the  problem  at  the 
end  of  the  first  iteration.  The  iteration  is  terminated  when 
the  nonlinear  solution  (frequencies)  obtained  during  two 
consecutive  iterations  differ  by  some  small  number  (say, 
one  percent). 

In  the  present  study  linear  (4-node)  and  quadratic  (8- 
node  and  9-node)  isoparametric  elements  are  employed.  The 
element  stiffness  matrices  are  of  the  order  20  x  20,  40  x  40, 
and  45  x  45  respectively. 


As  pointed  out  in  [30,61],  the  shear  deformable  theory 
presented  herein  can  be  derived  from  the  classical  thin  plate 
theory  by  using  the  penalty  function  method  to  incorporate  the 
slope-displacement  relations 


as  constraints  into  the  variational  formulation  of  the  thin 
plate  theory.  H  is  well-known  that  reduced  integration  must 
be  used  in  the  numerical  solution  of  penalty  function 
problems.  In  other  words,  the  shear  energy  terms  (which 
correspond  to  the  penalty  functional)  in  the  element  matrices 


M, 
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The  stresses  were  computed  at  the  Gaussian  points  (see  Fig. 
1): 


w 


w(f.f.0)E3h3  X  10^ 


Oj^(A,A,  ±t/2) , 


(4)2  ^  Oy(A.A,±t/4).  T,y  =  (4)2  T,y(B,B,±t/2)  (18) 

h 


"y  -  '-a'  "y' 

^  1  *-  V  =  a§^  ^yzM) 

in  layer  2  where  A  and  B  are  the  Gauass-point  coordinates  (as 
a  fraction  of  sides  a  and  b)  given  below: 

2x2  linear  4x4  linear  2x2  quadratic  4x4  quadratic 


A 

B 


0.125 

0.375 


0.0625 

0.4375 


0.05283 

0.4472 


0.02642 

0.4736 


u*^l'v»0 


edgeCJk^ 

y 

,^EDGE® 

f  '  1  •  ' 

>  a  f  a  1 

QUARTER 

1 

K"*— -a — •■--1 

PLATE 

• 

1 

1  0  •  0  < 

1 

_ ; 

1 

v-<ky-0 

SYMMETRY  BOUNDARY 
CONDITIONS 


9-  4'NOOE  ELEMENT 
-•-B-NODE  ELEMENT 
•>o-9-N00E  ELEMENT 


Figure  1.  Laminated  plate  geometry,  boundary  conditions, 
and  finite  element  discretization  by  linear 
and  quadratic  elements. 


Table  1  Effect  of  reduced  integration  on  the  maximum  deflec¬ 
tion  and  stresses  of  a  three-ply  (0*/90*/0')  square 
plate  subjected  to  sinusoidal  loading  (hj=h3=h/4, 
h2=h/2,  material  I). 


a/h 

Source 

w 

\y 

^xz 

3-OES 

[ni 

7.434 

0.599 

0.403 

0.0276 

0.301 

0.196 

CFS 

:3i. 

6.627 

0.496 

0.359 

0.0240 

0.414 

0.128 

Ref. 

[27] 

6.299 

0.532 

0.307 

0.0250 

— 

— 

4Q8-R 

6.627 

0.495 

0.359 

0.0240 

0.414 

0.128 

4L-R 

6.599 

0.467 

0.347 

0.0227 

0.395 

0.208 

4L-F 

6.427 

0.451 

0.328 

0.0219 

347.5 

247.2 

10 

2Q8-R 

6.615 

0.484 

0.351 

0.0234 

0.404 

0.126 

2Q8-F 

6.605 

0.483 

0.349 

0.0234 

0.404 

0.126 

2Q9-R 

6.561 

0.480 

0.341 

0.0231 

0.400 

0.126 

2Q9-F 

6.551 

0.479 

0.340 

0.0231 

0.399 

0.126 

2Q9-FR 

6.558 

0.480 

0.341 

0.0230 

0.400 

0.126 

2L-R 

6.508 

0.380 

0.284 

0.0187 

0.355 

0.107 

2L-F 

5.901 

0.334 

0.245 

0.0163 

261.6 

261.4 

TT] 

5.173 

0.0230 

0.328 

0.156 

CFS 

[311 

4.911 

0.524 

0.294 

0.0219 

0.434 

0.108 

Ref. 

.27. 

4.847 

0.557 

0.307 

0.0231 

4.911 

0.524 

0.294 

0.0219 

0.434 

0.108 

4L-R 

4.863 

0.494 

0.278 

0.0207 

0.415 

0.103 

4L-F 

4.346 

0.437 

0.245 

0.0183 

115.5 

115.3 

20 

2Q8-R 

4.091 

0.511 

0.287 

0.0214 

0.424 

0.106 

2q8-F 

4.876 

0.508 

0.284 

0.0214 

0.424 

0.106 

2Q9-R 

4.847 

0.505 

0.278 

0.0213 

0.418 

0.106 

2Q9-F 

4.827 

0.502 

0.275 

0.0211 

0.418 

0.107 

2Q9-FR 

4.847 

0.505 

0.278 

0.0212 

0.417 

0.106 

2L-R 

4.712 

0.404 

0.229 

0.0170 

0.353 

0.089 

2L-F 

3.236 

0.265 

0.149 

0.0111 

70.55 

70.4 

.11. 

■  U.T/ 1  ■ 

■tIICKPa 

urn?" 

CFS 

.31. 

4.337 

0.535 

0.269 

0.0212 

0.442 

0.100 

Ref. 

:27: 

4.363 

0.566 

0.284 

0.0223 

— 

— 

4.336 

0.535 

0.269 

0.0212 

0.442 

0.100 

4L-R 

4.281 

0.505 

0.254 

0.0200 

0.422 

0.096 

100 

4L-F 

1.034 

0.120 

0.060 

0.0048 

9.327 

9.404 

2Q8-R 

4.319 

0.521 

0.262 

0.0206 

0.435 

0.102 

2Q8-F 

4.143 

0.490 

0.244 

0.0199 

0.430 

0.100 

2Q9-R 

4.301 

0.418 

0.255 

0.0207 

0.428 

0.099 

2Q9-F 

4.158 

0.490 

0.242 

0.0201 

0.427 

0.100 

2Q9-FR 

4.299 

0.518 

0.255 

0.0206 

0.426 

0.099 

2L-R 

4.107 

0.413 

0.208 

0.0164 

0.356 

0.082 

2L-F 

0.315 

0.003 

0.015 

0.0012 

2.173 

2.192 

classical  thin- 

plate  solution 

4.350 

0.539 

0.269 

0.213 

0.339 

0.138 
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In  Table  1,  R  denotes  reduced  inteqration,  F  =  full 
integration,  and  FR  =  full  integration  for  bending  terms,  and 
reduced  integration  for  the  shear  terms.  The  following 
observations  can  be  made  from  the  results  of  Table  1. 

1.  The  nine-node  element  gives  virtually  the  same 
results  for  full  (3x3  Gauss  rule),  reduced  (2x2  Gauss 
rule),  and  mixed  (3x3  and  2x2  Gauss  rules) 
integrations.  However,  the  results  attained  by  using 
the  reduced  integration  are  the  closest  of  all  three 
integrations  to  the  closed-form  solution  [31]. 

2.  The  nine-node  element  and  the  eight-node  element 
(with  reduced  integration)  give  almost  identical 
results  (contrary  to  the  belief  that  the  nine-node 
element  is  superior  to  the  eight-node  element)  for 
all  side-to-thickness  ratios. 

3.  Integration  rule  has  a  more  profound  effect  on  the 
accuracy  in  the  eight-node  element  than  in  the  nine- 
node  element.  Full  integration  gives  less  accurate 
results  than  the  reduced  integration,  and  the  error 
increases  with  the  side-to-thickness  ratio.  This 
implies  that  the  reduced  integration  is  a  must  for 
thin  plates,  as  generally  recognized  for  all  penalty- 
function  based  finite  element  models. 

4.  Full  integration  results  in,  relatively,  smaller 
errors  for  quadratic  elements  and  for  refined  meshes 
than  for  linear  elements  and/or  for  coarse  meshes. 
That  is,  the  error  between  the  solutions  obtained  by 
2Q8-F  and  2Q8-R  is  smaller  than  those  obtained  by  4L- 
F  and  4L-R,  and  the  error  between  the  solutions 
obtained  by  4L-F  and  4L-R  is  smaller  than  those 
obtained  by  2L-F  and  2L-R. 

5.  Numerical  convergence  of  the  element  is  clear  from 
the  results. 

It  should  be  pointed  out  that  the  present  finite-element 
solutions  are  not  expected  to  agree  exactly  with  the  3-D 
elasticity  solution,  because  the  element  is  based  on  the  YNS 
theory  rather  than  the  3-D  elasticity  equations.  The  accuracy 
of  the  YNS  theory  is  evident  from  the  results.  While  the  YNS 
theory  seems  to  predict  the  maximum  deflection  very  close  to 
that  given  by  the  3-0  elasticity  solution,  the  stresses, 
especially  the  transverse  shear  stresses,  are  not  in  good 
agreement  for  thick  plates.  The  stresses  predicted  by  the  YNS 
theory  are  on  the  lower  side  of  the  3-D  elasticity  results. 

In  the  results  to  be  presented,  reduced  inteqration  was 
used  with  linear  and  quadratic  elements. 


Table  2  Comparison  of  maximum  deflection  and  stressesof  a  three-ply  (0'/90°/0°)  square  plate  subjected 
to  sinusoidal  loadinq  (h^  =  h/3,  material  I). 


U 


*The  first  row  of  stresses  is  for  4x4  Gauss  points  and  the  second  row  is  for  2x2  Gauss  points. 
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In  Table  2,  the  3-D  elasticity  solution  of  Pagano  [10] 
the  closed-form  solution  [31]  and  4Q8  and  2Q8  finite- 
element  results  (FEM)  are  compared  for  three-layer 
(0"/90"/0°),  equal  thickness,  square  plate  (material  I),  The 
present  finite-element  results  are  in  excel  agreement  with  the 
closed-form  solution  of  the  equations  governing  the  YNS 
theory.  The  stresses  in  3-D  elasticity  solution  were  reported 
to  be  maximum  at  the  following  points: 

°x  ""  ,  Oy  =  ±h/4)  ,  T^y  =  Tj(y(a,b,.^) 

"'xz  "  ’  Vz  ^  Ty2(-|.b,0)  (20) 

Table  3  compares  the  maximum  deflection  and  stresses  in  a 
three-layer  (h^  =  hj  =  O.lh,  h2  =  0.8h),  square  sandwich  plate 
(0°/90'’/0‘’ )  subjected  to  sinusoidal  loading.  The  material  of 

Table  3  Comparison  of  maximum  deflection  and  stresses  in  a 
three-layer  sandwich  (0‘’/90°/0°)  plate  subjected  to 
sinusoidal  loading  (h,=h3=0,  h2=0.8h,  material  I, 
2Q8-R). 


a/h 

source 

w 

5x 

■^xy 

X 

Ty2 

3-DES  [10] 

3.278 

0.4517 

0.2403 

0.185 

0.140 

2 

CFS 

9.552 

0.718 

0.3622 

0.1954 

0.876 

0.252 

FEM 

9.543 

0.716 

0.3615 

0.1961 

0.877 

0.253 

3-DES  [10] 

1.556 

0.2595 

0.1437 

0.239 

0.107 

4 

CFS 

4.7672 

0.867 

0.1520 

0.0877 

0.993 

0.174 

FEM 

4.7606 

0.865 

0.1517 

0.8781 

0.994 

0.174 

3-DES  [10] 

« 

1.153 

0.1104 

0.0707 

0.300 

0.053 

10 

CFS 

1.560 

1.017 

0.0776 

0.0533 

1.110 

0.095 

FEM 

1.558 

1.015 

0.0774 

0.0535 

1.112 

0.095 

3-DES  [10] 

• 

1.110 

0.0700 

0.0511 

0.317 

0.036 

20 

CFS 

1.0524 

1.053 

0.0595 

0.0450 

1.140 

0.076 

FEM 

1.0503 

1.051 

0.0594 

0.0452 

1.140 

0.076 

3-DES  [10] 

1.064 

0.0541 

0.0425 

0.323 

0.031 

50 

CFS 

0.9064 

1.065 

0.0539 

0.0424 

1.148 

0.070 

FEM 

0.9040 

1.062 

0.0538 

0.0425 

1.151 

0.071 

3-OES  [10] 

1.098 

0.0550 

0.0437 

0.324 

0.029 

100 

CFS 

0.8852 

1.067 

0.0531 

0.0420 

1.149 

0.069 

FEM 

0.8822 

1 .063 

0.0530 

0.0421 

1.158 

0.072 

CPT 

- 

1.097 

0.0534 

0.0433 

0.324 

0.029 

the  face  sheets  (i.e.,  layers  1  and  3)  is  the  same  as  material 
I,  and  the  core  material  is  transversely  isotropic  with 
respect  to  z  and  is  characterized  by  the  following  properties 
(see  Pagano  [10]), 

E,/Ep  =1.0  ,  E:,/E,  =  12.5  (not  used  in  YNS  theory) 

(21) 

0  I 3/G2 3^ 1 • 5E 2  ,  G32"0.4  E2  ,  2~ ^ i 3~ ^2 

The  finite-element  results  are  in  good  agreement  with  the 
closed-form  solution;  however,  the  YNS  theory  seems  to  predict 
the  stresses  quite  low  for  thick  plates. 

Figure  2  shows  the  nondimensional ized  deflection  versus 
side-to-thickness  ratio  for  two,  four  and  sixteen  layer  angle- 
ply  (45°/-45‘‘/+/-/ . . . )  square  plates  (material  II)  under 
sinusoidal  loading.  The  following  edge  conditions  were 
emp 1 oyed ; 


a/h  10  20  30  40  50 


Figure  2.  Nondimensional Ized  deflection  versus  side- 
to-thickness  ratio  for  simply-supported 
(SS-2)  angle-ply  plates  under  sinusoidal 
loading  (material  II). 
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Uo(o>y)  = 

=  U(j(a,y) 

=  0  , 

N6(0.y) 

=  0 

Vo(x,0)  = 

=  Vo(x,b) 

=  0  , 

N6(x,0)  ^ 

=  N^(x,b) 

=  0 

w(x,0)  = 

w(x,b)  = 

w(0,y) 

=  w(a,y) 

=  0 

(22) 

i|<x(x,0)  = 

=  0  , 

M2(x,0)  ^ 

=  M^(x,b) 

=  0 

1'y(0,y)  = 

=  'l'y(a,y) 

=  0  , 

Mi{0,y)  ^ 

=  0 

Dark  lines  and  symbols  correspond,  respectively,  to  the 
closed-form  solution  and  the  finite-element  solution  (with  the 
symmetry  conditions  of  the  CFS  imposed).  The  finite  element 
solutions  are  in  excellent  agreement  with  the  closed  form 
solution. 


Having  established  the  credibility  of  the  finite  element 
presented  herein  for  the  linear  analysis  of  layered  composite 
plates,  we  now  employ  the  element  in  the  nonlinear  analyses. 
First,  results  are  presented  for  single-layer  isotropic  square 
plates  under  uniform  loading.  The  essential  boundary 
conditions  used  are  {BC3  and  BC5  in  Table  4): 

simply-supported  (SS-3):  u=v=w=0  on  all  edges.  (22) 


clamped  (CC-1): 


u=v=w=0  on  all  edges, 

i|)y=0  along  edges 

parallel  to  x-axis,  (23) 

4»x=0  along  edges 
parallel  to  y-axis. 


Figures  3  and  4  show  the  nondimensional ized  center  deflection, 
w  =  w/h,  and  nondimensional ized  center  stress,  5  =  aa^/Eh^,  as 
a  function  of  the  load  parameter,  P  =  Pga'^/Eh'*  for  simply- 
supported  (SS-3)  square  plate,  under  uniformly  distributed 
load.  Figure  5  shows  similar  results  for  clamped  (CC-1) 
square  plate  under  uniformly  distributed  load.  The  results 
are  compared  with  the  Ritz  solution  of  Way  [65],  double 
Fourier-ser ies  solution  of  Levy  [66],  the  finite-difference 
solution  of  Wang  [67],  the  Galerkin  solution  of  Yamaki  [68], 
and  the  displacement  finite-element  solution  of  Kawai  and 
Yoshimura  [69].  Finite-element  solutions  were  computed  for 
the  five  degrees  of  freedom  (NOF  =  5),  and  for  three  degrees 
of  freedom  (NDF  =3);  in  the  latter  case,  the  in-plane 
displacements  were  suppressed.  Since  suppressing  the  in-plane 
displacements  stiffens  the  plate,  the  deflections  are  smaller 
and  stresses  are  larger  than  those  obtained  by  including  the 
in-plane  displacements.  Solutions  of  the  other  investigators 
were  read  from  the  graphs  presented  in  their  papers.  The 
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Figure 


Figure 


Load-deflection  curves  by  various  investigators 
for  simply-supported  (SS-3)  square  plate  under 
uniformly  distributed  load. 


p  so  100  ISO  zoo  zso 

P  •  q^(a/h)^/E 
0 


.  Nondimensionalized  center  stress  versus  the 
load  parameter  for  simply  supported  (SS-3) 
square  plate  under  uniformly  distributed 
load. 
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Figure  5.  Nondimensional ized  center  deflection  and  stress 
versus  the  load  parameter  for  clamped  (CC-1) 
square  plate  under  uniformly  distributed  load. 
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present  solutions  are  in  good  agreement  with  the  results  of 
other  investigators. 


Table  4 

Types  of  boundary  conditions 
study. 

used  in  the  present 

Notation 

(Type) 

Description  of  essential 

boundary  conditions. 

side  1 

side  2 

BCl 

(SS-1) 

V  =  w  =  lily  =  0 

u  =  w  =  \|»j^  =  0 

BC2 

(SS-2) 

U  =  W  =  tjly  =  0 

V  =  w  =  ijij^  =  0 

BC3 

(SS-3) 

u  =  V  =  w  =  0 

u  =  V  =  w  =  0 

BC4 

(SS-4) 

u  =  v  =  w  =iiiy  =  0 

u  =  v  =  w=  =  0 

BC5 

(CC-l) 

u  =  v  =  w  =  i|i^  =  0 

u  =  v  =  w=iiiy  =  0 

BC6 

(CC-2) 

U  =  V=  W=i|;j^  =  i|iy  =  0 

u  =  v  =  w=  iii^  =  ii<y  =  0 

BC7 

(CC-3) 

u  =  w  =  i|ij^  =  0 

V  =  W  =  liy  =  0 

BC8 

(CC-4) 

w  =  =  0 

W  =  Illy  =  0 

Table  5  shows  nondimensional ized  center  deflection  w  = 
w  /h,  center  stress  5"  =  o^a^/Eh^,  and  edge  stress  5°  = 
Oja^/Eh^  ofaclamped  tcC-2)  square  plate  under  uniformly 
distributed  load,  q_.  The  boundary  conditions  are  of  type 
BC6, 


CC-2;  u  =  v  =  w=\|»^=<|»y  =  0  along  all  edges.  (24) 

The  present  solution  is  obtained  using  a  2  by  2  uniform  mesh 
(in  quarter  plate)  of  nine-node  isoparametric  elements  (2Q9) 
with  (R)  and  without  (F)  reduced  integration.  The  present 
solution  for  center  deflection  and  stresses  agree  very  closely 
with  the  finite  element  solution  of  Pica  et  al.  [70],  and  the 
analytical  solution  of  Levy  [71].  The  edge  stress,  for  some 
reason,  does  not  agree  with  the  other  two  results. 


Table  5  Nondimensional ized  center  deflection  (is),  center  stress  (o^),  and  edqe  stress  (o®)  for  clamped 
(CC-2)  square  plate  under  uniformly  distributed  load  (Qq). 
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Next,  results  of  larqe-def lection  bending  of  layered 
composite  plates  are  presented.  First  a  two-layer  (-e°/o°) 
angle-ply  square  plate  (material  II)  under  uniformly 
distributed  load  and  simply  supported  boundary  conditions 
(BC2)  was  analyzed.  Figure  6  shows  a  plot  of 
nondimensional ized  center  deflection  (w_/h)  and  stress  a  = 
Oj^a^/Ejh^)  versus  the  orientation  { 6)  of  the  layers  for  the 
load  parameter  value  (P  =  qQa'+/Ejh‘*)  of  625.  It  is 
interesting  to  note  that  the  linear  deflection  and  the 
nonlinear  deflection  curves  do  not  resemble  each  other.  More 
specifically,  the  linear  deflection  is  maximum  at  0  =  15° 
whereas  the  nonlinear  deflection  is  maximum  at  6  =  45°.  The 
results  were  obtained  using  a  2  by  2  mesh  of  nine-node  element 
with  reduced  integration  (2Q9). 


Figure  6.  Nondimensional ized  center  deflection  and  stress 
versus  the  orientation  0  for  twr-layer  angle-ply 
square  plate  under  uniformly  distributed  loading 
(material  II,  a/h  =  100,  BC2,  P  =  625). 

Figure  7  shows  nondimensional ized  deflection  (obtained 
using  2Q9  mesh)  versus  the  load  parameter  P  for  a  single-layer 
(0°)  orthotropic  (material  II)  square  plate  under  two  loadings 
(UDL  and  SSL)  and  two  boundary  conditions  (BC2  and  BC6).  The 
large-deflection  results  of  Chia  and  Prabhakara  [41]  cannot  be 
compared  with  the  finite  element  solution  because  the  boundary 
conditions  used  there  involve  specifying  the  in-plane  stress 
resultants  (Np  N2  and  N^).  Figure  8  shows  the  load- 
deflection  curves  for  four-layer  cross-ply  and  angle-ply 
square  plates  (material  II,  BC6)  under  uniformly  distributed 
load.  The  figure  also  contains,  in  the  case  of  the  angle-ply 
plate,  results  for  a  thick  (a/h  =  10)  plate.  For  the  load 
parameter  value  P  =  3500,  the  thick-plate  deflection  is  about 
25%  higher  than  the  thin-plate  deflection;  this  shows  that  the 
shear  deformation  has  significant  effect  on  the  deflections 
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Figure  7.  Load-deflection  curves  for  single-layer  (0®) 
orthotropic  (material  II)  plates. 


computed.  The  difference  between  the  thin-  and  thick-plate 
deflections  increases  with  the  load  parameter. 

Figure  9  shows  the  non-dimensional ized  deflection  versus 
the  load  parameter  for  two-,  and  six-layer,  anti -symmetric 
(0*/90”/0”/. . . )  cross-ply  rectangular  plates  (material  II) 
subjected  to  uniform  loading.  The  plate  is  assumed  to  be 
clamped  (CC-3),  as  described  in  BC7  of  Table  4.  The  present 
solution  (obtained  using  2Q8  mesh)  is  in  good  agreement,  for 
various  aspect  ratios,  with  the  perturbation  solution  of  Chia 
and  Prabhakara  [41].  Due  to  lack  of  tabulated  results  in 
[41],  the  relative  differences  in  the  two  solutions  cannot  be 
discussed. 

Figure  10  shows  similar  results  for  two-,  and  six-layer, 
angle-ply  (45*/-45"/-/+. . .),  clamped  (CC-2)  rectangular  plates 
(material  II)  subjected  to  uniform  loading.  Again,  the 
present  result  is  in  close  agreement  with  that  of  Chia  and 
Prabhakara  [41].  The  nondimensionalized  stress,  o^,  for  the 
cross-ply  and  angle-ply  plates  discussed  above  is  plotted 
against  the  load  parameter  in  Figure  11. 


The  effect  of  the  transverse  shear  strain  on  the 
deflection  and  stresses  on  the  load-deflection,  and  load- 
stress  curves  is  shown  in  Fig.  12.  Note  that  the  deflection 
for  a/h  =  10  is  about  30%  larger  than  that  for  a/h  =  100,  at  P 
=  10.  That  is,  the  deflections  predicted  by  the  classical 
thin-plate  theory  are  lower  than  those  predicted  by  the  shear 
deformable  theory. 


Figure  8.  Load-deflection  curves  for  four-layer  angle- 
ply  (45°/-45  745  7-45°)  and  cross-ply  (07 
90°/0°/90°)  square  plates  under  uniform  loading 
(material  II,  BC6). 


6.2  Free  vibrations  results  [30-32,44,60]. 

The  effect  of  reduced  integration,  and  the  eight-node 
element  versus  the  nine-node  element  were  examined  for  linear 
frequencies  of  natural  vibration,  using  a  three-layer  cross- 
ply  (0*/90*/0*)  square  plate  (material  II,  BCl).  Table  6 
shows  the  nondimensionalized  fundamental  natural  frequency  as 
a  function  of  side-to-thickness  ratio,  type  of  integration, 
and  type  of  element.  From  the  results  obtained,  it  is  clear 
that  the  reduced  integration  or  mixed  (or  selective) 
integration  yields  relatively  better  results  for  all  ratios  of 
side-to-thickness  ratio. 
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Figure  9.  Load-deflection  curves  for  two-layer  and  six- 
layer  cross-ply  (0°/90V. . . )  rectangular  plates 
(material  II,  BC7)  under  uniform  loading. 
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Figure  11.  Nondimensionalized  center  stress  versus  the  load 
parameter  for  (a)  cross-ply,  and  (b)  angle-ply 
plates  (material  II,  BC6). 
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Figure  12.  Nondimensionalized  center  deflection  and  stress 
versus  the  load  parameter  for  four-layer  (O^/M®/ 
gOVO")  simply  supported  (BCl)  square  plate 
(material  I)  under  uniform  loading. 


8lo«T'— tna^mrotnro 
LDrooorOf— <— 
i^a^ocoir>\or^oococo 


c^iCMCOO^o^or^^s.r^»— 
OOir)ro»r>f— ororo^ 
tr>^ooro»— cMf^i— oo»— 

Lno^ooJto^^^cocoa^ 


^f^oootn*— roooo 
a%<ocotf>a>r^coooLn 
I  roooroi— »— \oovDoo 

I  . . 

o^oc^Jtn^^^cococo 


^oroou>»-“ioooro»r> 
oouf>‘^ir>oi^ooa>tr> 
tn^oocof— cvj^oo'^co 

ir>a>oc\jtr>vor^cooo<y» 


CVJOOO*— 0»— rooro 

tfj^cooo*— Csjr>.i— oocsj 
ir)0>ocMir>\o»^ooooa> 


Kfr^r^ooo*— ^c\ja> 
oa>mroir>a^r^ooa>vD 
U^COOOOOf—f—  voovooo 

ir)^ocsjvn^r^oor^oo 


roo^'~vo»— r^a>oo^ 
0»— vo^lr>o^^coo^r^ 
in^oocoi— cj'^ovooo 

u^a>ocMtf)vor^oooooo 


cvj^msoocsjomoo 
^  csj  cvj  m  o 


Figure  13  shows  a  plot  of  the  nondimens i^ona  1  i zed 
fundamental  linear  frequency,  u  =  (ija^/h) /P/E2,  versus  the 
side-to-thickness  ratio  (a/h)  for  four-layer  cross-ply 
(0"/90"/90"/0'’),  and  eight-layer  angle-ply  (45°/45‘/+/- . . . ) 
square  plates  (material  II).  The  boundary  conditions  used  for 
the  cross-ply  plate  are  those  in  BCl  and  the  boundary 
conditions  used  for  the  angle-ply  plate  are  those  in  BC2.  The 
finite  element  solutions  are  compared  with  the  exact  closed- 
form  solutions  (see  Reddy  and  Chao  [31]).  The  finite  element 
solutions  were  obtained  using  a  2  by  2  mesh  of  eight-node 
quadratic  elements  (2Q8).  The  finite  element  results  agree 
very  closely  with  the  closed-form  solutions  in  both  cases. 
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Figure  i3.  Nondimensional ized  fundamental  linear  frequency 

versus  side-to-thickness  ratio  for  cross-ply  (BCl) 
and  angle-ply  (BC2)  rectangular  plates  (material  II) 


To  illustrate  the  accuracy  of  the  present  element  for 
nonlinear  vibration  of  plates,  linear  and  nonlinear 
fundamental  frequencies  obtained  by  the  present  element  are 
compared  first  with  those  available  in  the  literature.  Table 
7  contains  ratios  of  linear  to  nonlinear  fundamental 
frequencies  (<*>l/‘^l)  i'or  various  ampl itude-to-thickness  ratios 
(wg/h)  of  simply-supported  and  clamped  square  plates  (v  =  0.3) 
as  determined  by  various  investigators.  Since  the  other 
investigators'  analyses  were  based  on  von  Karman  plate  theory, 
not  including  the  in-plane  displacements,  we  set  b/h  =  1000 
and  suppressed  the  in-plane  degrees  of  frreedom  in  the  present 
model.  In  Ref.  [72]  a  mixed  (based  on  Reissner  type 
functional)  finite-element  model  was  employed,  while  in  Refs. 
[73]  and  [74]  the  cubic  displacement  (based  on  total  potential 
energy)  finite-element  model  was  used.  Chu  and  Herrmann  [75] 
and  Yamaki  [76]  employed  Galerkin  method  to  determine  the 
frequencies.  The  present  results  were  obtained  using  2x2  mesh 
of  nine-node  quadratic  element  in  the  quarter  plate,  whereas 


Table  7  Ratio  of  linear  to  nonlinear  fundamental 

frequencies  for  various  ampl itude-to-thick- 


ness  (c/h)  ratios  of  simply-supported  and 


clamped 

isotropic 

(v  =  0.3),  square 

plates 

frequency  ratio 

Source 

0.2 

0.4 

0.6 

0.8 

1.0 

IWI 

(2x2) 

0.9818 

0.8671 

0.7958 

0.7271 

172] 

(4x4) 

0.9818 

0.9331 

0.8670 

0.7957 

0.7270 

r73T~ 

(4x4) 

0.9818 

0.9331 

0.8670 

0.7958 

174] 

(5x5) 

0.9821 

0.9338 

0.8673 

0.7943 

0.7233 

75 

0.9809 

0.9297 

0.8702 

0.9853 

0.7131 

77 

iiwahh 

[eo] 

'2x2 

) 

0.9929 

0.9727 

0.9419 

0.9038 

0.8616 

r72r 

(4x4) 

0.9930 

0.9730 

0.9425 

0.9048 

0.8631 

1731“ 

(4x4) 

0.9930 

0.9731 

0.9427 

0.9052 

0.8637 

174] 

(5x5) 

0.9930 

TO 

0.9938 

0.9116 

0.8750 
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KEmm 

Table  8  Effect  of  side-to-thickness  ratio  (b/h)  and  edge 
conditions  on  nonlinear  frequencies  of  isotropic 
{v=0.3)  square  plates. 


Boundary  frequency  ratio 

Condition  XlxIO^  — o;2  U'4'  O  O  TTO” 


0.9499 

“OTO 

“03DS" 

0.8610 

orrgw 

^imo9 

EC22 

0.9638 

0.9818 

0.9331 

0.8671 

0.7958 

0.7271 

EC3 

0.9499 

0.9866 

0.9494 

0.8955 

0.8333 

0.7692 

EC4 

0.9638 

0.9871 

0.9512 

0.8990 

0.8383 

0.7754 

EC5 

0.9499 

0.9992 

0.9927 

0.9839 

0.9720 

0.9576 

ECl 

■aHi 

BB^ 

BBUia 

EC2 

2388.0 

0.9815 

0.9324 

0.8658 

0.7942 

0.7254 

EC3 

2332.3 

0.9862 

0.9478 

0.8925 

0.8290 

0.7639 

EC4 

2388.0 

0.9869 

0.9506 

0.8978 

0.8366 

0.7733 

EC5 

2332.3 

0.9992 

0.9927 

0.9838 

0.9720 

0.9576 

IAMBI 

BKiiliia 

Emuim 

EC2 

9308.7 

0.9808 

0.9300 

0.8620 

0.7894 

0.7203 

EC3 

8944.1 

0.9851 

0.9440 

0.8854 

0.8152 

0.8471 

EC4 

9308.7 

0.9864 

0.9487 

0.8943 

0.8315 

0.7669 

^convergence  criteria,  e  =  10"^ 

^these  aqree  with  those  of  Ref.  [70],  \i  =  wjjh/^ 


4x4  mesh  of  eiqht-node  quadratic  element  was  used  in  [72].  It 
is  clear  that  the  present  resutls  are  in  excellent  agreement 
with  the  other  approximate  solutions. 

Next  the  effect  of  boundary  conditions  and  thickness  on 
the  frequency  ratio,  is  investigated,  and  the  results 

are  shown  in  Table  8.  The  boundary  (or  edge)  conditions  used 
in  Table  8  are  listed  in  Table  4,  and  are  explained  below; 

ECl:  Three  (w, i()y) -degree  of  freedom  model,  with  w 

specified  to  be  zero  on  the  boundary  (i.e.,  simply- 
supported  edqes)-BC3(3D0F) 

EC2;  same  as  in  BCl,  except  the  tangential  rotations  are 
also  specified  to  be  zero  on  the  boundary  - 
BC4{3D0F) 

EC3;  Five  (u,v,w,  ij/y) -degree  of  freedom  model  with  u,v 
and  w  specified  to  be  zero  on  the  boundary  -  BC3. 

EC4;  Five-degree  of  freedom  model  with  the  transverse 
deflection,  in-plane  displacements  normal  to  the 
boundary,  and  tangential  rotations  are  specified  to 
be  zero  at  the  boundary  -  BC2. 

ECS:  Five-degree  of  freedom  model  with  only  the 

transverse  deflection  specified  to  be  zero  on  the 
boundary. 

As  can  be  seen  from  the  results,  nonlinear  frequencies 
are  more  sensitive  to  the  boundary  conditions  than  linear 
frequencies.  This  should  be  expected  since  the  coupling 
between  membrane  and  flexural  stiffnesses  is  not  present  in 
the  linear  theory.  Restraining  of  the  in-plane  displacements 
in  effect  stiffens  the  plate  and  conseuqently,  the  frequency 
of  vibration  increases.  Thus,  for  example,  the  (nonlinear) 
frequency  of  vibration  of  ECl  (in  which  u  and  v  are  zero 
everywhere)  is  larger  than  that  of  EC3,  although  the  linear 
frequencies  are  identical.  In  EC5  in-plane  displacements  are 
not  restrained  at  all,  and  therefore  thickness  has  no 
noticeable  effect  on  the  frequency  ratio;  EC4  is  a  five-degree 
of  freedom  version  of  EC2. 

Figure  14  shows  the  plot  of  frequency  ratio  versus  the 
ampl itude-to-thickness  ratio  for  various  boundary  conditions 
and  side-to-thickness  ratio  of  isotropic  (v  =  0.3)  square 
plates.  Since  the  plate  stiffness  is  proportional  to  the 
plate  thickness,  and  increases  with  edge  constraints,  the 
nonlinear  frequencies  are  greater  for  thick  plates  than  for 
thin  plates. 

Next,  large-deflection  free  vibration  of  orthotropic  and 
layered  composite  plates  is  discussed.  In  Table  9,  ratios  of 
linear  to  nonlinear  fundamental  frequencies  are  presented  for 
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Table  9  Linear  to  nonlinear  fundamental  frequency  ratio  for 
various  values  of  the  amplitude-to-thickness  ratio 
for  single-layer,  and  two-layer  cross-ply  square 
plates  under  various  boundary  conditions  (material: 
graph ite-epoxyj  Al  =  (Ji/p7F^). 


Boundary 

Condition 

amplitude-to-thickness 
0.2  0.4  0.6 

ratio  (w 
0.7 

o/h) 

1.0 

Linear 

frequency 

EC4 

(SL)1 

0.9792 

0.9227 

0.8432 

0.7555 

0.6713 

18.819 

EC4 

(TL-CP)2 

0.9801 

0.9262 

0.8510 

0.7681 

0.6865 

16.096 

ECS 

(TL-CP) 

0.9990 

0.9912 

0.9805 

0.9663 

0.9490 

11.122 

EC6 

(TL-CP) 

0.9990 

0.9912 

0.9806 

0.9664 

0.9491 

11.173 

EC4 

(SL) 

0.9643 

0.8736 

0.7611 

0.6523 

0.5581 

14.118 

EC4 

(TL-CP) 

0.9775 

0.9174 

0.8354 

0.7478 

0.6649 

13.956 

EC6 

(TL-CP) 

0.9976 

0.9907 

0.9797 

0.9651 

0.9478 

10.300 

^SL  =  Single-layer 

^TL-CP  =  Two-Layer,  Cross-Ply 


Figure  14.  Ratio  of  linear  to  nonlinear  fundamental  frequencies 
(y  ®  versus  the  amplitude-to-thickness  ratio 

(w  /h)^fof'‘‘isotropic  (v  =  0.3)  square  plates  (in-plane 
degrees  of  freedom  not  included). 
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Figure  15.  Effect  of  orthotropy  (6,2/E2  =  0.3846,  v,2  =  0.3) 
on  the  linear  nondimensionafized  fundamental  fre¬ 
quency  (A)  and  on  the  ratio  (w  =  )  of  linear 

to  nonlinear  fundamental  frequencies  or  square 
plates  (b/h  =  10). 


single-layer  orthotropic  and  two-layer  cross-ply  square  plates 
under  various  boundary  conditions.  The  material  properties 
used  are  those  of  graphite-epoxy  material  (Gj^  =  6^3  =  G^,), 
material  II,  and  the  edge  conditions  EC6  are  described  below: 

EC6:  same  as  EC4,  except  no  in-plane  displacements  are 
specified  on  the  boundary  (or  same  as  ECS,  except 
the  tangential  rotations  are  also  specified  on  the 
boundary) . 

Couple  of  observations  can  be  made  from  the  results  presented 
in  Table  9.  First,  there  is  no  difference  in  the  frequencies 
obtained  by  ECS  and  ECS,  In  other  words,  specification  of  the 
tangential  rotations  does  not  alter  the  vibration  frequencies 
appreciably.  As  one  might  suspect,  the  frequencies  obtained 
by  EC4  and  EC5  are  noticeably  different;  the  nonlinear 
fundamental  frequencies  of  EC4  are  much  larger  than  those  of 
ECS  and  EC6.  Second,  there  is  a  reduction  of  about  38%  in 
frequencies  from  single-layer  to  two-layer  (for  the  same  total 
thickness)  plates. 

Figure  15  shows  the  effect  of  orthotropy  (E1/E2  for  fixed 
"  0-3846,  Vj2  =  0.3)  on  the  linear  fundamental 
frequencies  and  on  the  ratio  of  linear  to  nonlinear 
frequencies  of  single-layer  and  two-layer  cross-ply  square 
plates  (b/h  =  10)  subjected  to  boundary  condition  EC4.  While 


(b)  Angie-ply  places  (45* /-45*) 


Figure  16.  Ratio  of  linear  to  nonlinear  fundamental  frequencies 
versus  amplitude  to  thickness  ratio  of  single-layer 
and  two-layer  square  plates  (E-|/E2  "  ^12/^2  " 

v^2  =  0.25). 
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the  linear  frequencies  increase,  the  ratio  of  linear  to 
nonlinear  frequencies  decreases  with  increasing  Ej/Eo. 

Although  the  linear  fundamental  frequency  of  the  angle-ply 
plate  is  larger  than  the  cross-ply  plate,  the  ratio 
smaller  for  angle-ply  plate  than  for  the  cross-ply  plate. 

This  indicates  that  the  nonlinearity  is  more  pronounced  in  the 
angle-ply  plate  than  in  the  cross-ply  plate. 

Lastly,  Figure  16  shows  plots  of  linear  to  nonlinear 
frequencies  versus  the  ampl itude-to-thickness  ratio  for  two- 
layer  clamped  (transverse  deflection  and  the  normal  rotation 
are  set  to  zero  on  the  boundary-BC8)  square  plates  of 
graphite-epoxy  material;  plots  for  thin  (b/h  =  1000)  and  thick 
(b/h  =  10),  cross-ply  as  well  as  angle-ply  plates  are  shown  in 
Fig.  16. 

7.  SUMMARY  AND  CONCLUSIONS 

A  finite  element  based  on  a  shear  deformable  theory  (see 
[14])  and  von  Karman  (large  rotations)  theory  of  layered 
composite  plates  is  reviewed  and  numerical  results  for  static 
bending  and  free  vibrations  are  presented  for  isotropic, 
orthotropic,  and  layered  composite  plates  of  rectangular 
geometry.  It  is  noted  that  the  nonlinearity  is  less 
pronounced  in  the  layered  composite  plates  than  in  the 
isotropic  plates. 

From  the  study  it  can  be  concluded  that  the  four-node  and 
nine-node  isoparametric  elements  (with  reduced  integration  for 
shear  energy  terms)  give  accurate  results  while  possessing 
simplicity  over  traditional  plate  elements.  As  pointed  out 
earlier  (also,  see  [79-81]),  the  element  can  be  viewed  as  one 
based  on  a  penalty  function  formulation  of  classical  thin 
plate  theory  with  the  slope  continuity  conditions  as 
constraints.  The  penalty  terms  (i.e.,  shear  energy  terms)  in 
the  stiffness  matrix  must  be  evaluated  using  reduced 
integration  in  order  to  avoid  so-called  locking  (i.e.,  to 
avoid  excessively  stiff  elements;  see  [82]).  An  extension  of 
the  penalty-function  method  to  enforce  the  inter-element  slope 
continuity  needs  further  investigation  (see  Haugeneder  [83]). 

A  study  of  the  mathematical  properties  (e.g.,  stability,  error 
estimates,  etc.)  of  the  plate  bending  element  presented  herein 
is  not  complete  (see,  Ohtake,  Oden,  and  Kikuchi  [84,85]  for 
some  preliminary  results). 

Although  the  shear  deformable  theory  of  layered  composite 
plates  presented  here  yields  acceptable  solutions  for  global 
response  of  the  plate,  locally  (i.e.,  point-wise)  the  theory 
does  not  predict  stresses  accurately.  The  questions  relating 
to  interlaminar  stresses,  edge  effects,  and  delamination  in 
composites  (see  [86-95])  can  be  addressed  only  when  a  higher 
order,  three-dimensional  theory  is  employed  (see  [96-100]). 

The  comments  made  above  also  apply  to  the  case  of  laminated 
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composite  shells.  Thus,  developments  in  computational 
mechanics  related  to  finite-element  analysis  of  plates  and 
shells  in  the  next  decade  would  largely  concern  with  the 
development  of  computationally  simple  elements  that  are 
capable  of  representing  accurately  physical  features  of  the 
phenomena  involved. 
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APPENDIX  A:  NOMENCLATURE 


Ai j,Bij,D^j 

a,b 

E1.E2 


c  c  c 

^12*^13*^23 


extensional,  f lexural-extensional  and 
flexural  stiffnesses  (i,j  =  1,2,6). 

plate  planform  dimensions  in  x,y 
directions,  respectively 

layer  elastic  moduli  in  directions 
along  fibers  and  normal  to  them, 
respectively 

boundary  segments  of  domain  R 

layer  in-plane  and  thickness  shear 
modul i 


h,hi 


total  thickness  of  the  plate;  thickness 
of  i-th  layer 


I 


rotatory  inertia  coefficient  per  unit 
midplane  area  of  lamina 


^i 


shear  correction  coefficients 
associated  with  the  yz  and  xz  planes, 
respectively  (i  =  1,2) 

element  stiffness  coefficients  (a,e  = 
1,2, . .  .,5;  i,j  =  1,2, .. . ,n) 


L 


total  number  of  layers  in  the  plate 


Mi,Ni 


stress  couple  and  stress  resultant, 
respectively  (i  =  1,2,6) 


IJ 


element  mass  coefficients  (3, B  = 
1,2, ...  ,5;  i,j  =  1,2, ... ,n) 


n 

P 


number  of  nodes  per  element 

laminate  normal  inertia  coefficient  per 
unit  midplane  area 

shear  stress  resultant  (i  =  1,2) 

plane-stress  reduced  stiffness 
coefficients  (i,j  =  1,2,6) 
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% 


intensity  of  transversely  distributed 
load 


R 


R;  Rp 


1J  ’  1J 


t 


laminate  rotatory-normal  couplinq 
inertia  coefficient  per  unit  midplane 
area 

domain  (midplane  of  the  plate);  typical 
element 

element  matrices  in  FEM  formulation 
(C,n,^  =  0,  x,y;  i,j  =  1,2, ...,n) 

time 


Ui,U2,U3 


displacement  components  in  x,y,z 
directions,  respectively 


u,v 


in-plane  displacements  in  x,y 
directions 


U,U^,  etc. 


typical  variable  and  its  nodal  value, 
respectively 


x,y,z 


position  coordinates  in  cartesian 
system 


{A} 

^i 

®m 

p('") 


column  jf  vector  of  generalized  nodal 
displacements 

strain  components  (i  =  1,2,..., 6) 

orientation  of  m-th  layer  (m  = 

1,2, ...,L) 

density  of  m-th  layer  (m  =  1,2,...,L) 


A 

''i 

<>i 


'J'x’'1'y 


nondimensional  fundamental  frequency 

stress  components  (i  =  1,2,. ..,6) 

finite-element  interpolation  functions 
(i  =  1,2, ...,n) 

bending  slope  (rotation)  functions 

plate  curvatures  (i  =  1,2,6) 

fundamental  frequency  of  free  vibration 
(oil  =  linear;  =  nonlinear) 


(1) 
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APPENDIX  B:  ELEMENTS  OF  STIFFNESS  AND  MASS  MATRICES 


stiffness  matrix: 


[K]  = 


t[k21]  t[k22] 


a2[k23] 


u[k25] 

At[K"']  AtCK^^]  ^,-^33^  ^[-^343  ^,-^35^ 

^j^3|-^3]  +pA[Kf]  +uA[K^^] 


,31- 


u[K^^] 

u[k2^] 

u[Kf] 


t[K^^]  tCK""]  A[K;["] 


42. 


.43. 


.51 


t[K^']  A[Kp 


+a2[k43] 

.53. 


wCK^'^] 


w[K^^] 


w[K^^] 


+a2[k53] 


mass  matrix: 


[M]  = 


P[S] 

[0] 

[0] 

R[s] 

[03 

[0] 

P[S] 

[0] 

[0] 

R[S3 

[0] 

[0] 

P[S] 

[0] 

[03 

R[s] 

[0] 

[0] 

I[S] 

[03 

[0] 

R[S] 

[0] 

[01 

I[S] 

aft 

The  matrix  coefficients  K.^  are  given  by 

*  J 


[K^^]  =  A^iCS^^]  .  Aig([S^y]  .  [S^yf )  ^  AggCS^^]. 

[K^^]  =  +  A2g[Syy3  +  AggCs''^]^  = 

[K^^]  =  a^iCrJ'']  +  ^  ^  ^ 

"  "  Agg([Rjy]T  +  [R];y])  =  i  [k31]T. 

[K^"^]  =  B^^[S'‘’‘]  +  B^g([S'‘y]+[S''^]^)  +  BggLS^y]  = 

[K^®]  =  +  B2g[Syy]  +  =  [K^l]^ 

[K^^]  =  ^22^S^^1  +  A2g([S''^]^  +  [S^^^])  +  AggCs’^'']. 
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[k23]  =  A^2'^Rjy]^A22[Rj^]  +  A2g{[Rj^]^  +  [rJ^]  +  [rJ]) 

+  B2g[Syy]  +  B^gCs'^'^]  +  EggCS^^^] 


[K^^] 


[Kf] 

[Kf] 


B22[S^^]  +  B2g([s''y]  +  [S^y]^)  + 

BggCs'^'']  =  [K^2]T 

AssCs^^]  >  A^gCcs^y]  +  [s^y]T)  +  A^^csyy], 


3((i.  34). 

^2 


»  A55[S^0]  +  A^gCSyO]  * 

[K^"^]  =  B^^CR^""]  +  B^g([Rx^]  +  +  [Ry’']) 

+  8^2^^^]^  +  B2g[Rj;y]  +  Bgg([Rf  ]  +  [rJ^]) 

=  2[Kf  ]^ 

[Kf  ]  =  A^gCS^®]  +  A^^[Sy°]  =  [Kf  f , 

[k|^]  =  ^  ^  B22CRf  ]  +  B2g([Rjy]^  + 

[RJ^]  +  [Rf  ])  +  Bgg([Rjy]^  +  [Rf  ]) 

[k'*'‘]  =  D„[s’''‘]  +  D,g[s''y]  +  [s'^yf )  +  DggEsyy] 

AggCS], 


+ 
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A  po^fion  of  ^*1  report  will  appear  as  an  ASME  paper  to  be  presented  at  the 
1981  ASME  Vibration  Conference,  in  Hartford,  Connecticut  during  September  15- 
20,  1981,  The  paper  will  appear  in  its  entirity  in  "Recent  Advances  in  Non¬ 
linear  Computational  Mechanics,"  Hinton  et  al .  (eds.).  Pine  Ridge  Press,  Swansea, 
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angle-ply,  composite  materials,  cross-ply,  fiber-reinforced  materials,  finite 
element  method,  free  vibration,  laminates,  large  deflections,  natural  fre¬ 
quencies,  nonlinear  frequencies,  plate  bending,  rectangular  plates,  rotatory 
inertia,  transverse  shear  deformation. 

M.  ABSTRACT  (CantHm*  an  rarmaa  ataa  II  nMMAAir  arE  lEaNMIp  Ap  AtoaA  mmiaan 

'^A  finite-element  analysis  is  presented  for  large-deflection  bending  and  large- 
amplitude  free  vibration  of  rectangular  plates  of  cross-ply  and  angle-ply  con¬ 
struction  under  various  edge  conditions  and  loadings.  The  finite  element 
employed  in  the  present  study  accounts  for  transverse  sheairand  rotatory  inert  a 
(in  the  Reissner-Mindlin  sense),  and  large  rotations  (in  the  von  Karman  sense). 
The  finite  element  solutions  are  found  to  agree  very  closely  with  results 
available  in  the  literature.  — _  ■ 
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